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The short-distance behavior of exclusive amphtudes and in particular of the hadronic 
form factors at large momentum transfer has attracted continuous interest in the past decade. 
The theoretical basis is provided by the factorization theorem within a convolution formalism 
The nonperturbative information on the quark structure of the hadron is contained in the 
distribution amplitude which is the integral over the transverse momenta of the 

hadron valence-quark wave function. In a physical gauge (e.g. = 0), it is the probability 
amplitude for the hadron to consist of valence quarks with longitudinal momentum fractions 
< Xi < 1, = 1 (in a poo frame) moving collinearly up to a scale Q^. Contributions 

from higher Fock states are suppressed by powers of the momentum transfer. 

Theoretical constraints on the shape of the hadron distribution amplitudes can be ob- 
tained from QCD sum rules p and lattice gauge theory . For the nucleon, various model 
distribution amplitudes have been proposed and detailed studies have been carried 

out More recently, the quark structure of the A"''(1232) and higher resonances 

has been investigated [jl2|,|13],|T5 and models for the valence part of their distribution ampli- 
tudes have been worked out []T3|,|IB . 



It was shown in |]r7| that the transition form factor G%j and the neutron magnetic form 
factor G^.j are anti-correlated so that one of those form factors tends to be large when the 
other is small. If one uses the CZ or the COZ [0 nucleon distribution amplitude, 
turns out to be nearly zero, whereas IG^I/G^^ < 0.5. This possibility would explain why 
the data on the ratio Fe2/ Fmi yield a value compatible with zero. On the other hand, the 
GS [H nucleon distribution amplitude gives by construction \G\.j\/G\j < 0.1 and a \G%j\ 
of the same order of magnitude as G^j. This behavior is consistent with the possibility 
that o"„, the electron-neutron differential cross section, is dominated by G^, while G^j is 
asymptotically small or equivalently that ^ at all values |TH]. These results 



were obtained under the assumption that the A distribution amplitude can be approximately 
represented by the symmetric part of the nucleon one. However, the anti-correlation pattern 
still holds also for more realistic models for the A distribution amplitude (see []T3| and below). 



From the experimental side 19], there is still no clear evidence whether Q levels 



off or drops to zero at large Q^. A recent analysis of the unpublished data of the SLAG 



experiment E133 points rather to the second possibility, though the precision of the data is 
yet not sufficient to single out one possibility over the other. 

The theoretical predictions for the normalization and the evolution of G\i are sen- 
sitive to the shape of the proton distribution amplitude. It is therefore essential to have 
reliable models for both the nucleon and the A distribution amplitudes to be able to make 
comparisons with the data in precise detail. 

In a previous letter we have shown that two widely used models for the nucleon dis- 
tribution amplitude, the Chernyak-Ogloblin-Zhitnitsky (COZ) model and the Gari-Stefanis 
(GS) model — treated previously as competing alternatives — can be unified into a single 
model: the heterotic model. This model provides the possibility to analyze the data on 
cr„ under the assumption that \G\i\I(j'm < 0.1 (like the GS model), while matching the 
sum-rule requirements on the moments of $Ar up to the third order with almost the same 
overall accuracy as with the GOZ model. Furthermore, the calculated decay widths of the 
charmonium states ^Si, ^Pi, and ^P^ into are in excellent agreement with the data pi 



without individual adjustment of the various parameters. No other existing model is so 
successful in that respect |TT|J2^J2^. 

In the present letter, we apply similar ideas to derive an optimum distribution amplitude 
for the A"*" (1232) isobar. The new element of our approach is that we treat the sum- 



rule analysis of Garlson and Poor in conjunction with that of Farrar et al. [|16|. A 
second result of this paper is the calculation of the transition form factor G\.i in remarkable 
agreement with the data. In addition, we make predictions for the decay widths of the 
charmonium states ^Si, '^Pi, and ^P2 into A A. 

To leading order, the quark distribution amplitude for the baryon with helicity +1/2 can 
be represented in the form: 

$[i/2](x„Q2) = $..(x.) f:p[i/2](/i2)f^4^y"i>„(x.) , (1) 

where ^asip^i) = 120xia;2a;3 and {$n(xj)} are the eigenfunctions of the interaction kernel 



of the evolution equation, expressed in terms of Appell polynomials [|T|. The correspond- 
ing eigenvalues 7„ equal the anomalous dimensions of the lowest-twist three-quark opera- 
tors carrying the appropriate baryonic quantum numbers and are perturbatively calculable 
renormalizat ion-group coefficients ||2^. Note that in the A case only eigenfunctions sym- 
metric under Xi ^ X3 contribute (i.e. = = 0). 

The expansion coefficients B^^"^^ are nonperturbative quantities representing matrix el- 
ements of three-quark operators interpolating between the vacuum and the helicity +1/2 
baryon (the nucleon (N) or the A^(1232) isobar), renormalized at a scale /x^: 

< n|O("i"2"^)(0)|6[^/2](p) /;V(A)(2 ■ J5)"i+"^+"^+^^t^/2l(p) 0("i"2n3)_ (2) 

Here z is a lightlike vector (2;^ = 0), with z ■ q = = q'^ + for any vector q, 5*1^/^' (p) 
is the spin function of the baryon with helicity +1/2, and /Ar(A) is a dimensionful constant 
denoting the value of the matrix element at the origin. 



To carry out sum-rule calculations P,^j6|JT0.15.16], one uses correlators involving two of 

the 0("i"2"3): 



j(nin2n3,m)(^^^) = t I (i^X 6^^'^ < ^]|^(0 ("1"^"^^ (0)6^^^ (x)) |^] > (^ " 7) 



77' 



= (2 . g) "■l+"'2+"3+ni+3j(nin2n3,m)^^2^ 

where [z ■ 7)77/ serves to project out the leading twist-structure of the correlator. To obtain 
constraints on the moments of the baryon distribution amplitude, 

^in.n^n,) ^ ^^^^ , (4) 

J 

{[dx] = dxidx2dx3 5(1 — xi — X2 — x^)) a short-distance operator product expansion is per- 
formed at some spacelike momentum /i^ where quark-hadron duality is valid. By virtue of 
the orthogonality of the eigenfunctions the coefficients i?jj^/^l(/i^) and the constant fN(A) 
can be determined by means of moments inversion using as constraints the sum-rule re- 
quirements. The baryon distribution amplitude is then expressed in the form of a truncated 
series of Appell polynomials of up to second order, which means that we take into account 
the first six terms: n = 0, 1, ... 5. 



The coefficients for the heterotic nucleoli distribution amphtude were determined in |P0|| : 
Bq = 1 (due to the normahzation condition /o^[(ix] ^^ixi) = I ), Bi = 3.4437, B2 = 1.5710, 
i?3 = 4.5937, -B4 = 29.3125, and B^ = —0.1250. [Here and below the conventions and 
analytical expressions given in |T0[ are used.] The profile of ^^eterotic jg gj^own in Fig. |1|. 

In this note we treat the sum rules defined in Eq. (3) for rii + ^2 + 7^3 < 3 and m = 1 
by exploiting the possibility of simultaneously satisfying the moment constraints on the A 
distribution amplitude of both the Carlson and Poor |]15| and the Farrar et al. [T^ analyses. 
In this way we obtain a A distribution amplitude which has the explicit form 



= {7.2041 - 8.2859x2 - 21.3682xiX3 - 4.9247(xf + x^)} (5) 

with coefficients B^ given in Table | in comparison with those we determined for the specific 
model forms proposed by the above authors. Remarkably, like the nucleon case, this solution 
has heterotic character (see Fig. H). 

From Table ^ we can see that all the moments of the heterotic solution are within the 
range calculated by FZOZ for the amplitude T/^{xi,X2,x^) = ^a{x2,X3,Xi). Furthermore, 
the CP constraints on the amplitude Va(xi, X2, X3) = ^a{xi,X2, X3) +Aa{xi, X2, X3) are also 
satisfied, with the exception of the moment V^^'^^^ for which the heterotic solution yields 
a value slightly smaller than the estimated minimum value of the corresponding sum rule. 
Note that if one takes the estimated margins of both analyses [^,0 as they stand, then it 
is not possible to find a solution for the A distribution amplitude of the form given by Eq. |l| 
which simultaneously satisfies all sum-rule constraints. 

Following [|12| we input the heterotic nucleon distribution amplitude to calculate the 



— A+ transition form factor G*^, modeling the A"*" isobar by the three options labeled 
CP, FZOZ, and Heterotic. The evolution is due to the one-loop approximation of OsiQ"^) 
with Aqcd = 180 MeV [^. Here the average of two coupling constants {Q"^) and 
afiQ"^) is used, with arguments weighted by the virtualities determined for each model. For 
the heterotic A distribution amplitude we obtain af'^^''^\Q'^) = [«s(Q^ x 0.494) ^^(Q^ x 
0.088)]^/^; the corresponding values of asiQ"^) for the other considered models are given in 



20| and references cited therein. 

The results are shown in Fig. ^, where also comparison is made with the predictions 
for G%i derived from <|)^eterotic conjunction with previous models for nucleon distribution 
amplitudes. The experimental data are compiled in [|19|. In all cases the CP value |/a| = 
11.5 X 10"'^ GeV^ has been used, which is within the spread of the FZOZ estimate. We 
emphasize that the sign of G^,^ predicted by CZ [|], COZ 0, and GS P comes out negative 
for all A distribution amplitudes discussed here (cf. ||22|). It is only the heterotic nucleon 
distribution amplitude and the KS p one that yield a positive sign for G\,j. 

In order to account for (unknown) confinement effects at low Q^, we saturate as by 
introducing an effective gluon mass: ds{Q^) ^ ctsiQ^ + 4mg(Q^)). The situation is illus- 
trated in Fig. 0. In contrast to other approaches of this type |]8|j26[|, we use a dynamical, i.e.. 



scale-dependent gluon mass derived by Cornwall p7 



Due to the positivity of the anomalous dimension of the mass operator, this gluon mass 
vanishes asymptotically. This soft behavior at short distances leaves the validity of the 
form-factor evolution at large momentum transfer virtually unaffected. In the fit shown in 
Fig. H (dashed-dotted line), we take rag = 380 MeV, which agrees with Cornwall's consistency 
relation nig/ Kqcd ~ 1-5 — 2.0. 

Referring to the same figure, we see that including the perturbative evolution of the 
coefficients (cf. Eq. (1)), it is sufficient to provide a good fit to the data above ^ 
3 GeV'^ /c^ (dashed line). At lower values, additional nonperturbative parameters have 
to be introduced in the way just described (e.g., effective parton masses, quark clustering 
etc.) to account for the limitations of the leading-order formalism. 

In conclusion of this work, we make predictions for the exclusive decays of the charmo- 
nium levels ^S*!, ^Pi, and into AA. We here follow |TT,5D|. The branching ratio of the 
decay of the J*"^ = 1^^ state into AA is given by 

BRi :^ " 1 ^ , \m^)\ (?) 



0.75 167r2 


/a 


ln(M/A) 729 





where M ^ 2mc ~ 3 GeV and A = 0.4 GeV (the last value is taken from ||2^-see also [| 
The analogous expression for the J^*" = 2++ state has the form 

~ n QKt \4 16 /a '"^ 

^ 0.85 TTds - ^ 

^ ^ 729 

state into A A is 



BRl 



A A' 



^ ^ a// 

The partial width of the J^'" = V 



A\2 



(M-) 



(8) 



r(=^^i ^ AA) 



TTO;., 



1280 1/^1' 


/a 


2437r M 





A\2 



(9) 



where determines the value of the ^S'l-state wave function at the origin. Its value can 



be extracted from the leptonic width Ti^Si 



e^e 



(4.72 ± 0.35) keV ill via the Van 



Royen-Weisskopf formula. The result is |/^| = 383 MeV with mj/^ equal to its experi- 
mental value The nonperturbative input is due to /a and the decay amplitudes M^, 
M^, and M^. The results for the specific models considered here are summarized in Ta- 
ble |T|. For the branching ratio of the ^Si state, the heterotic distribution amplitude yields 
BR{^Si AA/^Si all) = 0.30 x 10-2%, where the new g value Ttot = 85.5t|^ keV is 
used. Note that in all considered decays as(mc) = 0.210 ± 0.028 (see third paper of p8[] ). 

We have performed our analysis pretending that factorization applies to exclusive reac- 
tions, so that a perturbative treatment is justified at accessible momentum transfer. This 
issue has been questioned by Isgur and Llewellyn-Smith and also by Radyushkin []32 



who argued that soft contributions dominate even at enormous momentum scales, render- 
ing perturbative QCD inadequate for exclusive reactions. However, more recently, Li and 
Sterman |3^ have shown that there is a an infrared protection of the perturbative picture 
provided by Sudakov effect suppression. 

In summary, we believe that the heterotic distribution amplitudes ^^"^^ and ^^^^ are 
physically appropriate solutions in terms of which several exclusive reactions can be theo- 
retically described without invoking additional assumptions and without tuning the various 
physical parameters from case to case. 
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FIGURES 



FIG. 1. Profile of the heterotic nucleon distribution amplitude. 



FIG. 2. Profiles of the model distribution amplitudes for A+(1232) described in the text. 

FIG. 3. Comparison with available data of the transition form factor 7pA^ calculated with 
the heterotic nucleon distribution amplitude and three different model distribution amplitudes for 
the A"*" isobar, as explained in the text. The solid horizontal lines are predictions for the absolute 
value of G%j at Q"^ = IbGeV"^ /c^ of other nucleon distribution amplitudes labeled CZ Q, GS 
KS Q, and COZ The data denoted by open circles are from |14]. 



FIG. 4. Comparison with available data of the transition form factor 7pA"'" calculated with the 
heterotic nucleon distribution amplitude and the heterotic distribution amplitude for the A+ isobar. 
Three different calculations are shown: the one-loop approximation of as{Q'^) (solid line), a modi- 
fied expression for UsiQ"^) which takes into account a dynamical gluon mass nig 

(Q^) (dashed-dotted 

line) , and the effect of Q^-evolution of the coefficients (dashed line) . 
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TABLES 

TABLE I. Expansion coefficients, B^, of the model distribution ampfitudes for A+(1232), 

discussed in the text. 

CP FZOZ Heterotic 

B^ 0.35 -0.175 -0.2499 

Bi" 0.4095 1.071 0.3297 

5^ 0.1755 -0.486 -1.6205 
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TABLE II. Moments ni + 77-2 + ^3 < 3 of the model distribution amplitudes for A"*" (1232) 



in comparison 


with the sum-rule constraints 


1. The numbers 


in parentheses 


are those given by 


FZOZ 0. 










Moments 


Sum rules 




Models 




(nin2n3) 


Ta{FZOZ) 


CP 


FZOZ 


Heterotic 


(000) 


1 


1 


1 


1 


(100) 


0.31—0.35 


0.35 


0.325 (0.32) 


0.321 


(001) 


0.35—0.40 


0.30 


0.350 (0.36) 


0.357 


(200) 


0.14—0.16 


0.16 


0.150 


0.14 


(002) 


0.15—0.18 


0.123 


0.160 


0.151 


(110) 


0.07—0.1 


0.101 


0.080 (0.07) 


0.078 


(101) 


0.09—0.13 


0.089 


0.095 (0.1) 


0.103 


(300) 


0.06—0.09 


0.085 


0.083 (0.085) 


0.073 


(003) 


0.06—0.10 


0.060 


0.085 (0.081) 


0.071 


(210) 


0.025—0.04 


0.039 


0.030 (0.025) 


0.027 


(201) 


0.04—0.06 


0.035 


0.037 (0.04) 


0.040 


(102) 


0.035—0.06 


0.031 


0.037 (0.039) 


0.040 


Va{CP) 


(001) 


0.33—0.37 


0.35 


0.325 


0.321 


(002) 


0.14—0.18 


0.16 


0.15 


0.14 


(101) 


0.072—0.12 


0.095 


0.088 


0.091 
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TABLE III. Charmonium decays in AA for the models discussed in the text. 



Amphtudes CP FZOZ Heterotic 

11418.23 21585.99 11651.26 

Mi" 30924.07 48233.23 26277.40 

1480.67 1882.31 1134.94 

Observables 

0.311 X 10-3% 1.113 X 10-=^% 0.324 x 10-^% 

0.100 x 10-3% 0.106 X 10-3% 0.072 x 10-^% 

0.439 X 10-2 keV 0.709 x 10-^ keV 0.258 x 10"^ keV 



Vi^Si AA) 
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